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Abstract 

Interplay between kinetic roughening and phase ordering is studied in a growth SOS 
model with two kinds of particles and Ising-like interaction by Monte Carlo simu- 
lations. We found that, for a sufficiently large coupling, growth is strongly affected 
by interaction between species. Surface roughness increases rapidly with coupling. 
Scaling exponents for kinetic roughening are enhanced with respect to homoge- 
neous situation. Phase ordering which leads to the lamellar structure persisting for 
a long time is observed. Surface profiles in strong coupling regime have a saw-tooth 
form, with the correlation between the positions of local minima and the domain 
boundaries. 
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1 Introduction 



Growth by vapour deposition is a technologically important process for pro- 
ducing high quality materials. During last years much progress has been made 
in understanding of growth mechanisms on the microscopic level [1]. How- 
ever, these results refer mainly to single-component homogeneous growth and 

1 Electronic address: kotrla@fzu.cz 



Preprint submitted to Elsevier Preprint 



1 February 2008 



growth of binary materials is much less understood. Essential feature in bi- 
nary systems is phase ordering [2]. Recently there was interest in interplay 
between phase ordering during growth and kinetic roughening [3-5]. Kinetic 
roughening is one of the important aspects of crystal growth. Due to stochastic 
fluctuations during a deposition process the growing surface becomes rough 
and roughness typically increases with time as a power with a characteristic 
exponent (3. This is important from the practical point of view, besides kinetic 
roughening also attracted a lot of interest in the context of non-equilibrium 
statistical mechanics [6]. 

In this paper we discuss the interplay between compositional ordering and 
kinetic roughening in a simple model for growth of binary alloys which we 
recently formulated [5]. 



2 Model 

Our model is based on the single-step solid-on-solid geometry. We confine 
ourselves here to one-dimensional substrate with the coordinate i. The surface 
is described by a single- valued function hi(t) with the additional constraint 
\hi — hj\ = 1 for neighbouring sites % and j. A new particle can be added 
only at a growth site corresponding to a local minimum in the surface height 
profile. We consider two types of particles, which we distinguish by a variable 
a taking the value +1 or —1. Growth rules are controlled by a change of 
energy of alloy after deposition of a new particle. The energy is given by the 
Ising-like interaction. The probability of deposition of a new particle of type a 
to a growth site i is proportional to exp(—AE(i, a) / ksT); ks is Boltzmann's 
constant and T is temperature. The change of energy is AE(i,a)/ksT = 
—Kg [a(i — 1) + a(i) + a(i + 1)] — Ha. Here K is a dimensionless coupling 
strength and H is an external field. We call this model a two-component 
single-step (TCSS) model. We consider here only the case of zero external 



The surface roughness is described by the surface width w(t, L) = (y h 2 — h ), 
where t is the time, L is the linear size and the bar denotes a spatial average, 
(...) a statistical average. Ordering is studied by the time evolution of the 
average domain size along the surface D(t). 



3 Results 

We performed extensive simulations for various values of K. Evolution of the 
surface profile is affected by composition of the surface, in particular it de- 
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pends on composition of the initial flat substrate. We considered three different 
types: i) a neutral substrate, i.e. without any interaction with deposited par- 
ticles, ii) an alternating substrate, with alternating types of particles and iii) 
a homogeneous substrate composed of one type of particles. 

Evolution on the neutral and the alternating substrate is almost the same. In 
Fig. 1 we can see that the effective exponent {f3 e & = 0.45) is at first larger 
than the exponent for homogeneous growth, f3 = 1/3, which corresponds to 
the so-called Kardar-Parisi-Zhang (KPZ) class [7]. However, after a certain 
time t cr oss ~ 10 3 the effective exponent crosses over to the value close to 1/3 
(/3 e ff = 0.32). This behaviour can be explained by evolution of compositional 
ordering. In the inset of Fig. 1 we present time dependence of the average 
domain size for three substrate types considered. We can see that for all of 
them the domain size saturates at time close to t C ross- Hence, ordering of the 
surface leads to enhanced increase of the roughness, but when ordering stops 
the surface roughness increases further with the same exponent as in one- 
component single-step growth model. 

In the case of the homogeneous substrate the behaviour of the surface width 
changes dramatically (cf. Fig. 1). The general features are the following. At 
first the width starts to grow with the exponent close to 1/3 (this regime is 
extended over several decades of deposition time for very strong coupling). It is 
caused by the fact that practically no domains with minority type of particles 
are formed because probability to nucleate these domains is rather small; 
growth is the same as in the pure single-step model. Then the slope crosses over 
to a very large value (/? e ff = 1-15). This unexpected effect is related to growth 
of trapezoidal features with the length given by positions of the minority type 
domains, which start to appear but often again quickly disappear. The rapid 
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Fig. 1. Surface width w vs. time for coupling strength K = 1.0, and three different 
types of the substrate of size L = 10000. Data are averaged over 30 independent 
runs. The values of slope, (3 = 0.45 and (3 = 0.32, were obtained by fit to points for 
the alternating substrate, the straight line with slope 1/3 is only guide for eye. 
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increase of roughness ceases when relatively stable minority type domains 
are nucleated. The further increase of roughness is very slow up to the time 
when the average domain size saturates and growth continues with the KPZ 
exponent (3 = 1/3. We illustrate the difference in evolution on the neutral and 
on the homogeneous substrate in Fig. 2. We can see that for the homogeneous 
substrate dark domains appear very slowly. 

Let us now turn to the dependence on the coupling K in the case of the neutral 
substrate. For small K, there is no ordering in the growth direction. Small do- 
mains appear and shrink again, but for a larger K the clear lamellar structure 
exists. We illustrate this in Fig. 3 which shows evolution of configurations on 
a substrate with periodic domains imitating the lateral superlattice. We can 
see that for K = 2 lamellae persist for long time. Notice the correlations be- 
tween the positions of local minima of surface profile and domain boundaries. 
We observed that this is a typical case for stronger coupling independently 
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Fig. 2. Examples of evolution of surface profiles for K = 1.0, and two different 
substrates: (a) the neutral substrate and (b) the homogeneous substrate. Surface 
profiles at various times increasing as powers are indicated by black lines, only part 
of the grown material close to the surface is shown for given time. Dark and light 
grey correspond to different types of particles. System size is L = 250. 
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Fig. 3. Examples of evolution of surface profiles on the substrate with alternating 
domains of different particle types for several coupling strengths K. Meaning of 
colours and other conditions are the same as in Fig. 2. 

rapidly with the coupling. The effective exponent j3 e s for initial growth in- 
creases with K and it seems to approach the value 1/2; for K = 2 we have 
found (3 = 0.52 ± 0.02. The crossover time for change of /? eff to the KPZ 
exponent is also increasing rapidly with K (t cross oc e' yK with 7 pa 3.5). 



4 Conclusion 



We demonstrated that phase ordering leads to faster kinetic roughening than 
in the homogeneous case. However, after some time this behaviour crosses 
over to the ordinary behaviour for homogeneous growth. But the crossover 
time increases rapidly with interaction between particles, so practically only 
the regime with enhanced roughness can be seen for stronger coupling. We 
also illustrated the strong sensitivity of evolution on the initial composition 
of the substrate. 
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